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Inter arrival  Times 


Section  1. 

Introduction  and  Summary. 

A random  variable  X is  said  to  have  an  IMRL  ( increasing  mean 
residual  life)  distribution  on  [O^oo)  if  Pr(X  < O)  = 0^  Pr(X  > t)  > 0 
for  all  t,  EX  < and  E(x-t|x  > t)  is  increasing  in  t > 0*  A 
random  variable  X is  said  to  have  a DFR  (decreasing  failure  rate) 
distribution  on  [0^<»)  if  Pr(X  < O)  = 0^  Pr(X  > t)  > 0 for  all  t, 
and  the  conditional  distribution  of  X-t  given  X > t is  stochastically 
increasing  in  t > 0. 

Consider  a renewal  process  with  interarrival  time  distribution  F^ 
an  IMRL  distribution  on  [O^oo).  Let  M(t)  be  the  expected  number  of 
renewals  in  [O^t]  including  an  initial  renewal  at  time  zero.  If 
M2  = x^dP(x)  < 00  define  U(t)  = tf  We  show  (theorem  4) 

r 00  k+2  / \ 

that  if  ""  •'O  ^ dF(x)  < 00  for  an  integer  k > 0 then: 


(1) 


U(t)  > M(t)  > U(t) 


min 

0 < i < k 


c . t 
1 


In  (l)  c^  is  an  explicitly  computed  function  of  ‘ ‘ ’ ^±+2’ 

0^...^k.  Moreover  the  numbers  v.  = c /c . are  ipcreasing  and 


for  V.  T<t<Vo  c.t^=  min  c.t~^i  thus  for  te[v.  ^^v.] 
J-l-  - J 0 0<i<k  ^ -1'^  J 

the  lower  bound  in  (1)  reduces  to  U(t)-c ,t 

^ ^ a t 

We  also  show  (theorem  5)  that  if  ii''„(a  ) = / e ° dF(t)  < oo  for 

-bo  0 

an  a > 0 then  for  0 < a < a : 
o — o 


1 


(2)  U(t)  >M(t)  > U(t)-(e^^-l)■^[(^i^a)■^-(^J2/2^L^)-(ya)-l)"^] 

2 

where  as  before  U(t)  = t/ . 

The  bounds  (l)  and  (2)  give  intervals  for  M(t)  whose  lengths 
rapidly  approach  0 as  t 

We  also  derive  several  monotonicity  results  for  IMRL  and  DFR 
renewal  processes  (theorems  2 and  3)-  If  F is  IMRL  then  the  expected 
forward  reciirrence  time  EZ(t)  is  increasing  in  t as  is  M(t)-t/ii^^ 
and  M(t+h)-M(t)  > h/p.^  for  all  t>0,h>0.  If  F is  DFR  then 
A(t)  (the  renewal  age  at  time  t)  and  Z(t)  are  stochastically  increasing 
in  t>  N(t+h)-N(t)  is  stochastically  decreasing  in  t for  all  h > 0; 
M(t+h)-M(t)|.h/p^  as  t -» ooj  M(t)  is  concave  and  if  F is  absolutely 
continuous  then  the  renewal  density  decreases  to  as  t -»oo. 

Our  results  follow  from  a representation  theorem  (theorem  l)  for 
stationary  IMRL  renewal  processes,  and  general  delayed  DFR  renewal 
processes.  This  representation  seems  to  be  an  ideal  tool  for  the  study 
of  IMRL  and  DFR  renewal  processes. 

The  theory  of  DFR  distributions  is  developed  in  Barlow,  Marshall  and 
Proschan  [5],  Barlow  and  Marshall  ([2], [3])^  Barlow  [l],  and  Barlow  and 
Proschan  [4].  II/IRL  processes  are  studied  by  Bryson  and  Siddiqui  [6],  and 
Haines  and  Singpurwalla  [8]. 
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Section  2 . 


Definitions • 

A probability  distribution  F is  said  to  be  DFR  on  [O^oo)  if 

F(0  ) = 0,  F(t)  < 1 for  all  t,  and  F(t+s)/F(t)  is  increasing 

in  t > 0 for  all  s > 0 (F(x)=1-F(x) ) . By  increasing  (decreasing) 

we  mean  monotone  non-decreasing  (non-increasing).  If  F(0  ) 0, 

F(t)  < 1 for  all  t}  and  F is  absolutely  continuous  with  pdf  f, 

then  F ■ is  DFR  on  [0,oo)  if  and  only  if  there  exists  a version  of 

f for  which  h(x)  = f(x)/F(x)  is  decreasing  ([5]p*378).  The  function 

h is  called  the  hazard  function. 

If  jj.,  = / *”  xdF(x)  < “ then  define  G(x)  = 11,^  F(y)dy.  G 

■^0  0 

is  the  stationary  forward  recurrence  time  and  renewal  age  distribution 
for  a renewal  process  with  interarrival  time  distribution  F.  A prob- 
ability distribution  F is  said  to  be  IMRL  on  [O^m)  if  F(0  ) = 0^ 

F(t)  > 0 for  all  t,  < <»,  and  E(X-t|x  > t)  = p,^G(t)/F(t)  is 
increasing  in  t.  Since  F(t)/|i^G(t)  is  the  hazard  function  of  the 
distribution  G we  see  that  F IMRL  <=>  G DFR.  Also  note  that  if 
F is  DFR  and  |j,^  < « then  F is  IMRL.  It  is  easy  to  construct 
examples  for  which  F is  IMRL  but  not  DFR. 

Let  Xq=  0^  and  X^^X^^--.  be  i.i.d.  with  distribution  F.  Define 
n 

Sp,  = 0^  S = L X.^  n=1^2,...,  R(t)  = < t]  = min[i:S.  > t]  ^ and 

M(t)  = EW(t).  Define  F(t)  = l-F(t),  q = F(o),  n = /”  x^dF(x) 

2 WCt)  ° 

L(t)  = M(t)-t/ Zi(t)  = E X^-t  (the  forward  recurrence  time 

N(t)-1 

at  t)  and  A(t)=t-  E X.  (the  renewal  age  at  t).  The  process 

n 1 
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(N(t)j  t > 03  is  called  an  ordinary  F renewal  process.  If 

are  independent^  Pr(X!,  > t)  = /°°  (F(t+y)/F(y))dH(y)  where  H is  a 

^ 0 

probability  distribution  on  [0^oo)j  Xj  ~ F for  j > 2^  ~ ^ ^1'’ 

N'(t)  = < t]  = max{i:  S'^  < t} , then  {N'(t),  t > O]  is  called 

a delayed  F renewal  process  with  initial  age  distribution  H.  For  a 
delayed  F renewal  process  Z'(t)^  Z’(t),  N’(t)^  M'(t)^  X^^  will 
denote  the  anologues  of  Z(t)j  A(t),  W(t),  M(t),  X^>  S^.  A delayed  F 
renewal  process  with  initial  age  distribution  G(x)  = F(y)dy 

is  called  a stationary  F renewal  process.  A stationary  F renewal 
process  satisfies  EN'(t)  = t/p^  and  Z'(t)  ~ A'(t)  ~ G for  all 
t > 0. 
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Section  3 « 


Representations . 

We  will  either  assume  F IMRL  on  [0,co)  or  DPR  on  [0^«>).  We 
will  construct  two  dependent  renewal  processes.  Process  1 will  be  an 
ordinary  F renewal  process.  Process  2 will  be  a stationary  F renewal 
process  in  the  IMRL  case,  and  delayed  F renewal  process  with  arbitrary 
initial  age  distribution  in  the  DFR  case.  The  special  feature  of  process 
2,  under  this  construction,  is  that  S„,.  , = S'. , i=l,2,  ...  for  a random 
integer  W (the  distribution  of  W will  depend  on  the  initial  age 
distribution).  Processes  1 and  2 differ  only  in  that  process  2 has  zero 
renewals  in  [0,S^)  while  process  1 has  W renewals  in  this  interval. 

The  simple  nature  of  this  difference  is  exploited  to  obtain  our  results. 

The  construction  is  based  on  a simple  idea  which  is  obscured  by  the 
details  of  the  construction  and  proof.  We  decompose  the  hazard  into  2 
components, the  first  component  causing  failure  for  both  processes,  the 
second  component  only  causing  failure  for  process  2.  The  construction 
uses  the  following  lemma: 

Lemma  1.  Let  X be  distributed  as  F where  F will  either  be 

assumed  IMRL  on  [0,oo)  or  DFR  on  [0,co).  Set  K(t)  = G(t)  in 

the  IMRL  case  and  K(t)  = / °°  (f( t+y)/F(y) )dH(y ),  where  H is  an  arbitrary 

0 _ _ _ 
probability  distribution  on  [0,co),  in  the  DFR  case.  Define  " K(t+v)/K(v), 

J.y.('t)  = F(t)/K^(t).  Then  is  the  survival  function  of  a perhaps 

defective  distribution  on  [0,oo). 

Proof . In  the  IMRL  case  = (F(t)/G(t)  )(G(t)/G(t+v))G(v),  and  since 

G is  DFR  both  F(t)/G(t)  and  G(t)/G(t+v)  are  decreasing.  Thus 


j 


is  decreasing.  In  the  DFR  case  J (t)  = K( v)[ /'”  (F(t+v-ty)/l'^t)) (F(y) ) dH(y)'] 

^ 0 

Since  F is  DFR  the  denominator  is  increasing^  and  then  since  the 

numerator  is  constant,  is  decreasing.  Thus  in  both  cases 

is  decreasing.  In  addition  is  right  continuous,  equals  1 for 

t < 0 and  is  always  between  0 and  1.  It  is  thus  the  survivad 

function  of  a perhaps  defective  distribution  on  [0,<»).  j| 

We  proceed  with  the  construction.  Again  K(t)  = G(t)  when 

F is  assumed  IMRL  , and  K(t)  = /*”  (F(t+y)/F(y) )dH(y)  with  H 

0 

an  arbitrary  probability  distribution  on  [0,oo)^  when  F is  assumed  DFR. 


Construct  and  W^  independent  with  ~ K,  W^  ~ J where 

J(t)  = F(t)/K(t).  If  Z,  < W-,  set  X-,  = X!  = Z^  and  X.  = X'  = Y.  ,, 

' 1-1  111  j j j-r 

j=2,3^..*  where  > l)  is  an  i.i.d.  sequence  with  distribution 

F independent  of  (Z^,W^) . If  2,^  > W^,  set  X^  = W^  and  go  to  stage 
2 . At  stage  2 construct  Zg  and  W^  conditionally  independent  of  each 
other  and  of  (Z^,W^)  given  W^,  with  = v having  distribution 

K^(t)  = K(t+v)/K(v),  and  W^lw^  = V distribution  J^(t)  = F(t)/K^(t). 
If  Z„  < W_  then  set  X„  = Z„,  XI  = W^+Z^  and  X.  = X'  = Y. 

2—2  2 2112  j j-1  j-2 

j=5,4,...  where  {Y^,i=l,2, . . .}  is  i.i.d.  with  distribution  F and 
independent  of  (Z^,  W^,  Z^^  W^).  If  W^  < set  and  go  to 

stage  5-  We  reach  stage  m if  and  only  if  W^  < Z^,  i=l, ...,m-l,  in 


which  case  X.  = W.,  i=l, ...,m-l.  At  stage  m we  construct  Z and 
11  m 

W^  conditionally  independent  of  each  other  and  of  (Z^,W^)  • • '(Z^  l-’^m  1^ 
m-1  m-1  m-1 

given  Z W.,  with  (Z  | Z W.  = v)  ~ K , (W  | Z W.  = v)  ~ J • 


If  Z < W set  X = Z , X'  = Z W.+Z  , X.  = X'.  = Y.  , j=m+l,  . 

m—  m m m l ^^im  j j-m+1  j-m 

where  {Y^,i  > l]  is  i.i.d.  with  distribution  F and  independent  of 
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(Z^  ,W,  • ,Z  , ,W  If  Z >W  we  go  to  stage  m+l  and 

1-^  1.  2 2’  ’ m-l  m-l  m m ^ 

repeat • 

Theorem  1.  Under  the  above  construction: 

(i)  i=lj2,  • . . } are  i-i.d.  with  distribution  F. 

(ii)  [X'oi  > l)  are  independent^  X'  ~ X'.  ~ F for  j=2^3* 

1 -L  J 

(lii)  X'.  = X,,,.  T for  i=1^2^*»*  where  N = minti:  Z.  < W. } 

1 K+i-l  1—1 

and  Pr(lil  < oo)  = i. 


Proof,  (i)  Define  K = minli:  Z.  <W.},N=a>  if  W.  <Z.  for  all 
1 — 1 11 

i.  Now  (X  I N > i,  (¥.,Z  ) = (w  ),  J=l^...^i-l)  ~ mln(Z*^W*) 

J J J J V V 

i-1  ^ ^ * 

where  v=  Z w,Z  ~KjW  ~J  and  Z and  W are  independent. 

]_^"vv"vv  V V ^ 

Since  K^(t)j^(t)  = F(t),  (X.|n  > i,  (W  ,Z  ) = (w..z.)  j=l,  ~ F 


Moreover  (X. |N=j  < i)  ~ Y.  . ~ F.  Thus  X.  ~ F independent  of 

X 1 J 1 

X,j...^X._,  . Since  this  holds  for  all  i,  (X.^j  > 1}  is  i.i.d.  with 

X 1 X J 

distribution  F. 


(ii)  In  our  construction  we  generated  It  will  now  be 

convenient  to  continue  constructing  W.'s  for  j > N.  At  stage  j 

J 

construct  W.  to  be  conditionally  independent  of  W given 

J X J —x 

j -1  I J "1 

ZW.^  with  (W.  Z W.=v)  ~ J . We  know  that  Pr(Z  W.=oo)  = 1 because 
1 1 ^ j'  1 1 ^ V 11 

K is  DFR^  so  inf  Pr(W.  > t)  > inf  Pr{W  > t)  = F(t)/lim  K (t)  > 0. 

i V V-^oo 

Thus^  given  t,  for  almost  all  (W^=w^^  i=1^2^...}  we  can  find  j so 
J-1  j 

that  Z w.  < t < Z w,  . Then 
11  - 1 1 


j-1 


j-1 


Pr(X'  > t|w,,Wp;,  ...)  = ( n K (w.))  K , (t-  Z w ) = K(t)  . 

1 Id  n 1 (?-!  ^ J;;!  1 1 


£=1 


Z w^i 
1 ^ 


Z w. 
1 ^ 


Thus  X^  ~ K. 
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w 

(iii)  Since  X'  ~ 

1 1 

renewals  in  a finite  time 


and  process  1 can  have  only  finitely 
interval,  and  since  Pr(X^  < <»)  = 1,  it 


many 

follows 


that  Pr(N  < oo)  = 1.  | 
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Section  4. 


Some  properties  of  IMRL  and  DFR  renewal  processes . 

Theorems  2 and  3 below  present  several  properties  of  IMRL  and  DFR 
renewal  processes.  Theorem  2 is  extensively  used  in  section  5 to  obtain 
oiur  bomd  for  M(t). 

We  will  need  a simple  result  which  is  well  known^  but  for  which  we 
have  no  reference. 

Lemma  2.  Let  F be  a distribution  on  [O^co).  if  < oo  define 
F^(t)  = /”F(x)dxi  if  jjg  < 00  define  FgCt)  = / F^(x)dx.  Then: 

( i ) [j,^  ^ 00  ^ tVt  ) -.0  as  t -» <»>  for  k > 0^  *\+i  ^ “ 


t^(t)dt  = [Jjj+p/ ^■‘■1  t^^(t)  0 as  t oo;  for 

^>0,  < “ =>  /“'t^3^(t)dt  = < “ and 


as  t ->  ooj  for  k > 0^ 

^.J^^y(k+l)(k+2)(k+3)  < <». 


< “ ^ /r  t^^(t)dt  = 


(ii)  For  a > 0^  ^-nCa)  = / e^^dF(t)  < ■»  implies  e^^F(t)^ 

^ 0 

e F,  (t)  and  e FQ(t)  converge  to  0 as  t -^co.  Moreover 


/q  e^^F(t)dt  = a ■^(  \l;p(  a) -l)  < «>,  /“  e'^'^F^(t)dt  = a ‘^(  a)-ap^-l)  < 


00  at— 


and  /“  e^^F2(t)dt  = a a)-a^(  P2/2)-ap^-l)  < 


Proof.  The  equalities  between  integrals  follow  from  interchanging  the 
order  of  integration.  The  convergence  to  zero  follows  from  the  equality 
between  the  integrals  in  each  pair  and  integration  by  parts.  || 


Theorem  2 . Consider  an  ordinary  F renewal  process  with  F IMRL.  Then: 

(i)  M(t)-t/|i  and  EZ(t)  are  increasing  in  t > Oj  M(t+h)-M(t)  > h/p. 
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if  iJ,^  < “ 


for  all  t > 0,  h > 0 and  converges  to  h as  t -4  oo; 
then  L(t)  = M(t)-t/|j,^-lJ2/^|J.^  t 0 as  t -><». 

(ii)  If  iJg  “ "then  0 > L(t)  > / (G(x)-g  ^(x))dx 


> 


“1 

*^1  G(x)dx^  where  q = F(o). 


(iii)  For  k > Oj  lim  t^(t)=0  and  < “ 

00  It  “^00  ^ , 

=#>  0 > / t L(t)dt  > -CO  j For  a > 0,  ^„(a)  - J dF(t)  < oo 

0 ^0 
^ lim  e®'^L(t)  = 0 and  0 > / e^^L(t)dt  > -<»  . 


t -4  00 


0 


(iv)  If  ij^<  00  and  h is  measurable;,  boionded^  and  lim  h(t)  = 0 


tig- 

then  lim  [f^  h(t-x)dM(x)-ii^^  h(x)dx]  = 0. 


t — >00 


Proof . Recall  that  process  1 is  an  ordinary  F renewal  process^ 
process  2 a stationary  F renewal  process,  and  i=l>***  • 

Note  that  by  Wald’s  identity  EN  = whether  or  not  is 

finite  • 


(i)  M(t)-t/M^  = E(R(t)-N'(t)).  Since  N(t)-R’(t)tR 

M(t)-t/ix,  t EN  = 1x5/2 by  the  monotone  convergence  theorem.  Thus 
^ 4 N(t) 

if  1X5  < 00  then  0>L(t)l0  as  t->oo.  Since  Z(t)  = E X.-t  it 

1 ^ 

follows  from  Wald’ s identity  that  EZ(t)  = |x^(M(t)-t/ix^),  thus  EZ(t)t  • 
Since  M(t+h)-M( t )-h/p^  =--  (M( t+h)-( t+h)/p^)-(M( t)-t/p^)  and  M(x)-x/p^ 
is  increasing,  M(t+h)-M(t)  > h/p^.  The  convergence  of  M( t+h ) -M( t ) 
to  h/pj^  can  be  proved  di.rectly  from  the  construction,  but  we  will  simply 
appeal  to  Blackwell's  theorem  noting  that  if  F is  a lattice  distribution 
with  period  a>  then  E(X-oy2|x  > (a/ Z)  = E(X-ay'2|x  > O) 

= E(X|X  > 0)-<u|2  < E(x|x  >0)  so  F is  not  IMRL-. 
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(il)  By  the  argiiment  in  (i)  L(t)  = n^^(z(t)-Z' (t)) . By  construction 


0 > z(t)-z'  (t) 


4 


r 0 if  t > x[ 


Z(t)-(X[-t)  if  t < X^  . 


Thus  li^L(t)  = G(t)  E(Z(t)|x^  > t)  - /”  G(x)dx;.  Since  Z(t)|x^  > t 

is  a mixture  of  distributions  of  the  form  X-v|x  > v where  X ~ F 

-1 


which  is  IMRL,  E(z(t)|x|  > t)  > E(x|x  > O)  = q 


M-1’ 


Thus 


0 > L(t)  > J°°  (G(x)-q  V(x))dx  > /“  G(x)(3x  . 

ij 

(iii)  These  results  follow  from  the  inequality  0 > L(t)  > f°  G(x)dx 


aad  lecima  2* 


(iv)  Consider  h( t-x)dM(x)-u,^  h(x)dx  = /”  h(t-x)l(x  < t). 

0 -^0  0 ~ 


d(M(x)-x/|j^)^  where  l(x  < t) 


0 

1 if  X < t 
0 if  X > t 


By  assumption  the  integrand  converges  to  zero.  By  part  (i)  of  this 
theorem  M(x)-x/ii^  t Pg/Sp^  as  x thus  if  pg  < «>  M(x)-x/p^ 


2 -1 


Thus 


has  total  finite  variation  on  [O^oo)  equal  to  p^/Bp^-q 

h(t-x)l(x  < t)  is  dominated  by  sup|h(s)!,  an<i  I sup  h(  s ) | d(M(x)-x/pL.  ) 

2 s ■ ■ s 

= ( sup|h(  s)  I )(  p /2p, -q  ) < 00.  The  result  follows  from  the  dominated 


convergence  theorem. 


Theorem  3 • Let  F be  DFR  on  [O^oo).  Then; 

(i)  A(t)  and  Z(t)  are  stochastically  increasing  in  t. 
( ii  ) N(t+h)-N(t)  is  stochastically  decreasing  in  t and 

M(t+h)-M(  t)  ^ L/p^  as  t ->oo. 

(iii)  M(t)  is  concave. 
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(iv)  If  F is  absolutely  continuous  then  as  t ->  oo^ 

where  m is  the  renewal  density  function. 

Proof,  (i)  Give  process  2 initial  age  distribution  H = F„  where  F„ 
is  the  age  distribution  of  process  1 at  time  s.  Then  since  A*(t)  > A(t) 
and  Z'(t)  > Z(t)  for  all  t by  construction,  Pr(A(t+s)  > a = Pr(A'(t)  > a) 

> Pr(A(t)  > a)  and  similarly  for  Z(t). 

(ii)  Pr(N(t+h)-N(t ) > k)  = /^q  * Pr(N(h-x)  > k).  Since 

Z(t)  is  stochastically  increasing  and  Pr(N(h~x)  > k)  is  decreasing  in 
X,  Pr(N(t+h)-N(t ) > k)  is  decreasing  in  t.  Since  N(t+h)-N(t)  is 
stochastically  decreasing,  M(t+h)-M( t)  is  decreasing. 

(iii)  We  want  to  show  that  for  0<x<y,  0<o;<l, 

M(o:x+(l-Q:)y)  >0:  M(x)+( l-0;)M(y),  equivalently  that  r(x,(l-Q:)(y-x)) 

> r(x+( l-o:)(y-x),  0!(y-x))  where  r(t,s)  =■  (m(  t+s)-M(t))/s . Since  M 

is  right  continuous  (M  is  actually  continuous)  it  will  suffice  to  show 
that  r(tj^,k^/n)  > r(t2,k2/n)  for  all  n,  0 < t^^  < t^  + k^/n  < tg. 
consider  the  case  0 < k^  < kg,  thfe  case  k^  > kg  > 0 follows  similarly. 

By  part  (i)  of  this  theorem  r(tg,k^/n)  < r(tj^,k-|/n)  and 
r(tg+k^/n,(kg-ki)/n)  < r(tg+k^/n,l/n)  < r(t^+(k^-l)/n,l/n)  < r(tj^,k^/n), 
thus  r(tg,kg/n)  = ( k^/kg )r( t2,k^/n)+[ ( kg-k^)/kg]r( tg+k^/n, ( kg-k^)/n) 

< (k^/kg)r(t^,k/n)+((k2-k^)/kg)r(t^,k/n)  = r(t^,k/n). 

(iv)  Let  h(x)  = f(x)/F(x)  the  decreasing  hazard  function.  Now 
m(t)  = Eh(A^)  with  h^  and  A(t)  stochastically  increasing.  Thus 
m(t)  is  decreasing  and  therefore  has  a limit j by  the  elementary  renewal 
theorem  this  limit  must  be 
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Section  5 • 


Bounds  for  M(t). 

Theorem.  If  F Is  IMRL  on  [O^oo)  and  V^+2  ^ °°  ^ integer  k > 0 

then: 


(3) 


U(t)  > M(t)  > U(t)  - min  c.t 

0<i<k  ^ 


2 2 

where  U(t)  = t •’ 

0 < = -i  s^  ^(s)ds  = / s^d(M(s)-s/|i^)  > -oo  for  1=1,  ...^k 

The  term  is  a function  of  |j,^  • • • i=l,  ...,k  which  can  he 

recursively  computed  from: 

(4)  c^  = 7^-li^^i'  X 

s=l 

where  i+l)(i+2)p^]-[  |j2ia^^^/2(i+l)  |j.^] 


/ s • ) X- . - , i — 1 , . • • , k 

3'  i+l-s 


Equation  (4)  can  be  explicitly  solved  yielding: 


£ 

where  = [{±,...,1^):  > 1;>  r=l,...,£,  i^=k)  . 

Proof.  Since  M(t)  = U(t)  + L(t)  and  L(t)  < 0 (theorem  2,  (i)), 
M(t)  < U(t).  Since  L(t)  is  increasing  (theorem  2,(i)), 

L(t)  > L(0)  = M(0)-|^/2ia^  = = -c^.  Thus  M(t)  > U(t)-9Q. 


13 


Since  L(t)  is  increasing  it  is  at  least  as  big  as  its  average  over  [O^t] 

with  respect  to  any  probability  measvire  on  [0,t]-  Thus  L(t)>it  ^ s^  ^(s)dE 

—i  i“l  / \ 

which  since  L < 0 exceeds  it  J s Lis/ds  = -c.  > oo  by  theorem  2 


part  (iii).  The  equivalence  between  -it  ^ / s^  ^L(s)ds  and 

0 

/ s^d(M(  s ) -s/ [i^)  follows  by  integration  by  parts  and  part  (iv)  of 
theorem  2. 

To  identify  the  c.'s^  start  with  the  identity  M(t)  = 1 + M(t-x)dP(x)^ 

subtract  U(t)  from  both  sldes^  multiply  both  sides  by  t^  and  use  the 

identity  t^-(t-x)^  - Y.  (^)('t“x)^x^  This  yields: 

r=0  ^ 


(6)  t^L(t)  = / (t-x)^L(t-x)dF(x)  + h(t) 

0 

where  h(t)  = h.,  (t)+hp(t)-h,(  t);  li-|(t)  = /^  [ Z (^)(t-x)^x^"^]  • 

I ‘d  1 Q r 

L(t-x)dF(x),  ^^1  F(x)dx^  lij(t)  = t ( iJ2/2ia^)F(  t)  . 

Wow  (6)  is  the  renewal  equation  g = h+g*F  = h*M^  with  g(t)  = t^L(t). 

By  part  (iv)  of  theorem  2,  if  we  can  show  that  h(t)  is  bounded^ 

integrable  and  that  lim  h(t)  = 0,  then  we  can  conclude  that 

t 00 

lim  t^L(t)  = / h(t)dt.  But  for  i=l^...^k^  lim  t^L(t)  = 0 

t ->  oo  0 t -^00 

by  part  (iii)  of  theorem  2-  Thus  the  conclusion  will  reduce  to 


00 

/q  h(t)dt  = 0^  which  will  provide  us  with  a useful  Identity. 

To  show  that  h is  bounded^  integrable  and  convergent  to  zero^  we 
do  so  separately  for  each  h_.  . By  lemma  2,  h^  and  h^  are  convergent 
to  zero  and  integrable  with: 


ih 


(7) 


/“  h2(t)dt  = 


(8) 


/”  h^(t)dt  = p,^_^^y.2/2(i+l)ii^ 


The  bomdedness  of  h^.  and  h^  follows  from  the  boundedness  on 

finite  intervals  and  the  convergence  to  0 as  t ->oo. 

Defining  L(y)  = 0 for  y < 0 we  write 
“ i-1 

hnCt)  - f ( S (^)x^  (t-x)^L(  t-x)  ))dF(x)  . Since  (t-x)^L(t-x)  ->  0 

1 0 I--0  ^ 

as  t 00  by  part  (iii)  of  theorem  2)  the  integrand  converges  point- 
wise  to  0.  Moreover^  since  s^|l(s)|  is  bounded  on  finite  intervals 

and  converges  to  0^  sup  s^|l(3)|  < oo^  thus  the  integrand  is  dominated 

s . 

by  the  integrable  function  Y ( )(sup  s |l(s)|)x  . Thus  by  the 

r^  ^ s 

dominated  convergence  theorem  h,  (t)  ^0  . The  above  argument  also  shows 

i-1  . 

that  |h^(t)|  < (^)p^_^(sup  s^|l(s)|)  < <»_,  thus  h^  is  bounded,  h^ 

r=0  ^ ^ ^ s . ^ 

is  integrable  by  part  (iii)  of  theorem  2 since  / jh(t)|dt  = 

. , 0 
X-1  . 00 

S (^)p.  / s|L(s)|ds<oo.  Moreover: 

r=0  I”  X'X"  Q 


(9) 


00  ^ ^ r 

/ h (t)dt  = X (^)p.  I s L(s)ds  = -i'.  X (c  /r!  )X. 
0 j,_Q  ^ 0 ^ 


r=l 


i+l-r 


OO  0®  « 

The  identity  h(t)dt  = 0 is  equivalent  to  - h^(t)dt  = Jq  h2(t)dt  ■ 

/“  h (t)dt.  Using  (7)^  (8)^  and  (9)  this  gives  us  i'.  S (c  /r.')>>- 
0 ^ !r=-l  ^ i‘ -L“r  X 1 

reduces  to  (k) . 


which 
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Define  d„  = 0^  d. 

0 1 

1=1:,  • • - Mq  0:,  = 


i=l>  . 
i=l^ 


•^k.  Then  rewrite  (4)  as: 


(10) 


d.  = 5. 
1 1 


E <1  p. 

$ ^x-e 

s=0 


(lO)  is  a discrete  renewal  equation.  Its  solution  is  d.  = S .M.  . 

1 -^-n  ,1 

“ ,(s)  , 


where  M.  = ^ P. 


(s) 


where  P' 


.(s) 


s=0 


j=o  J i-J 

is  the  s-th  convolution  of  p.  Since 
i i-1 


= Oj  = 0 for  s > i.  We  thus  obtain  d.  = E & .M.  . = 6.  + Y B .M. 


1 
i-1 


j=i  ^ 


j=i 


j 1- 


. - S ^ 1 (TT  ),  which  is  equivalent  to  (5)* 

^ A ^(i,...,iJeA,_,  , 1 


At  first  glance  it  may  seem  that  we  need  to  compute  c^t 


several  1 in  order  to  compute  inf  c.t 


-i 


for 

Fortunately  this  is  not 


the  case. 

Lemma  3 « Assume  that  F is  IMRL  on  [O^oo).  Then  either  F is  quasi- 
exponential  (F(t)  = qe  j0<q<l,  A,>0^t>0)  in  which  case 
c^  = 0 or  F is  not  quasi -exponential  in  which  case  c^  > 0 for  all 
Define  v_^  = 1^  v^  = x-0 ,1)2, . . . , with  o/ 0 = 0^  oo/oo 


Then  v.  is  increasing  and  for  v.  .,  < t < v.,  c .t  ^ = inf  c.t 
^ J"-*-  J J i-0^  1;  . . . ^ 


Proof.  If  U-,  = “ then  r d(M(t)-t/iu^)  = 00  ^ so  C^  = 00^  i=0;l;...  . 

If  pg  < 00  then  by  theorem  2 part  (ii)  L(t)  is  the  distribution  function 

of  a positive  measure.  Thus  c^  =.0  if  and  only  if  L is  a constant.  Since 

lim  L(t)  = 0,  that  constant  must  be  0.  Let  Z'(t)  denote  the  forward 
t — > 00 

recurrence  time  of  the  stationary  process  at  time  t.  By  our  construction 
Z'(t)  > Z(t)  with  equality  for  all  t with  probability  one  iff 
~ x|x  > 0.  Since  L(t)  = ij^^(Z(t )-Z' ( t ) ) we  see  that 
c^  = 0 <=J>  G(t)  = q ^(t)  <=^  h*(t)  = qp^^  where  h*  is  the  hazard 


i6 


function  for  the  distribution  G.  But,  under  this  last  condition 

G(t)  = e , so  F(t)  = qe  ^ , thus  F is  quasi -exponential. 

2 

Next  consider  the  Hilbert  space  L [[0,oo),p,L]  where  P is 
the  collection  of  Borel  sets  on  [0,oo).  Kecall  again  that  L is 
the  distribution  function  of  a positive  measure.  Now 

e] 

by  Schwartz's  inequality.  Note  that  the  inequality  trivially  holds 

in  the  quasi -exponential  case,  and  also  holds  in  the  alternative 

case  since  c > 0 and  c.  = cxd  =>  c = oo.  The  inequality  is 

equivalent  to  v = c /c.  > = v.  (again  o|0  = 0, 

oo|oo  = oo).  Thus  the  v^'s  are  increasing. 

We  next  show  that  t < v.  implies  c.t  ^ for 

-3  ^ 3 - D-hn 


m — 2 ^ • • 

Thus  for  t < V.,  c.t”^ 

- J 3 

with  k > j. 

A similar  argument  >iiich 

■"  "i  “k 

t > V.  , implies  c.t  < c t for 

- J -1  3 - k 

that  c.t  ^ 

3 

is  optimal  in  [ v.  ,v.]. 

3 3 

m > 1.  Then 

c.t"^  < c.  if 

3 - 3tm 

m-1 

t”  < c.^  /c.  = Tf"  v.^.  . 

Since  t < v.  and  each  of  the  m terms  of  the  product  are  equal  to 
3 

or  greater  than  v.  (since  t)  it  follows  that 

3 3 


< TT  v,,, 
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Remark.  The  above  proof  shows  that  if  we  only  consider  c^,...,c  , 
0 k 

“k  “k. 

where  k is  finite  then  c,  t = min  c.t  for  all  t > v,  . 

k „ 1 — k-1 

0 < 1 < k , 

- - a„t 

Theorem  5*  If  F is  IMRL  on  [O^oo)  and  \lr.„(a^)  = /“  e dF(t)  < oo 
■ I?  u u 


for  an  > 0,  then  for  0 < a < a^: 


(11)  U(t)  > M(t)  > U(t)-(e^^-l)  [(|i^a)“^-(p2/2p^)^(tp(a)-l)“^]  . 


Proof.  The  proof  is  very  similar  to  that  of  theorem  4.  Choose 
ae(0,aQ].  Using  L(t)  < 0,  L(t)t  as  in  the  proof  of  theorem  3 
we  obtain: 


(12)  U(t)  > M(t)  > U(t)  + a(e^’*^-l)“^  /“  e^®L(s)ds 


00  0^-^ 

where  0 > /q  e L(t)dt  > -00  by  part  (iii)  of  theorem  2. 

To  evaluate  ^^^(a)  = Jq  e L(s)ds  we  start  with 

M(t)  = 1 . /*  M(t-.)aF(K),  subtract  U(t)  from  each  siae  ana 
multiply  both  sides  by  e . This  gives: 


(13) 


e®'^L(t)  = e®'^^“^^L(t-x)dF'(x)  + I (t) 


where 


l(t)  = , .gj_(t)  = /^(e^-l)e®'^^“^^L(t-x)dF(x);, 


'^2*'^^  = 1^^^  / F(x)dx,  ij(t)  = -(p^2iJ.^)e®'^F(t) 
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We  verify  the  conditions  of  paxt  (iv)  of  theorem  2 in  a similar 

manner  as  in  the  proof  of  theorem  making  heavy  use  of  lemma  2 

and  part  (iii)  of  theorem  2.  The  conclusion  of  part  (iv)  of 

theorem  2.,  in  li^t  of  lim  e L(t)  = 0 (part  (iii)  of  theorem  2 

t 00 


becomes : 


(l4)  (t^(a)-l)\lr^(a)=(p2y^M-^)[  (tj.(a)-l)/a]-(p^a^)  (\Jf^(a)-ti^a-l)  . 


Since  \|r„(a)-l  ^ 0 for  a ^ 0 we  can  divide  both  sides  of 

r 

(l4)  by  y„(a)-l  and  solve  for  ilf^(a).  This  gives: 
r L 


(15: 


\|/^(a)  = ) +[a(tj,(a)-l)] 


-1 


Substituting  (I5)  into  (12 ) gives  us  (11).  || 

Remsirk.  (ll ) and  (15 ) will  hold  for  a < 0 whether  or  not 

\lrp(a)  < 00  for  an  a > 0.  If  |a^  < 00  and  we  let  a t 0 (H) 

then  we  obtain  M(t)  > U(t)-c^t  ^ where  c^  is  given  in  (4).  In 

general  “ implies  (o"")  exists  and  equals 

-(k+l)"^  '^1  . 1 • Thus  tT(a)  can  be  considered  a generating  function 

for  the  c.'s.  However  unless  the  particular  form  of  t„(a)  leads 

1 r 

to  a simple  expression  for  (a),  expressions  (4)  and  (5)  of 


theorem  3 will  be  preferable  for  computing  the  c^’s. 
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Section  6. 


Improved  bounds  'when  F is  DFR. 


The  bounds  given  in  theorem  k,  corollary  3}  and  theorem  5 for 
IMRL  distributions  can  be  improved  for  DFR  distributions.  Define 
«„  = 1,  a.  = (i/i+i)^.  i > 1,  c*  = a.  c.  -tdiere  c.  is  given  in 
theorem  3>  and  = c*,,/c*  = (a  /a.  )v. . Also  define 
g (t)  = a(e^^-l)~^  se^^ds  = (te^^/e^^-l)-a  and  i|f(a)  = -aijf  (a) 

a 0 b 

[ (l+j_a)~^-(p2/2|i^)-('ip(a)-l)“^]. 

Corollary  U.  Assume  that  F is  DFR  on  [O^oo).  Then: 

(i  j If  p.  “ then  U(t)  > M(t)  > U(t)-  min  c*t"^ 

“ “ 0<i<k  ^ 

(ii)  t and  for  v^  ^ < t < c*t"'^  = inf  c'!ft"^  : 

thus  for  V*  ..  < f < V*  the  bound  in  (i)  is  given  by  U(t)-c*t"^*. 

i “b  0 J 


(iii)  If  tT,(a^)  < 00  for  an  a^  > 0 then  for  0 < a < a^: 

F 0 0 — 0 

ag"^(t)  -1  , , -1 

U(t)  > M(t)  > U(t)-(e  -1)  ^(a)  > U(t)'*(e®'^  "^-l)  itf(a) 


Proof.  (i)  ii  is  concave  by  theorem  3 part  (iv).  Thus 
/q  s-s^-^ds)  = L((j/J+I)t)  > /q  s^“^L(s)ds 


> Iq  s^~h(s)ds  = Thus  L(t)  > -c^[(a+l/j)t]''^’ 

V — n 

= c.t  • The  argument  now  proceeds  as  in  theorem  3- 


(ii ) A simple  differentiation  argument  shows  that  • 

Since  by  corollary  5 and  v*_^^/vf  = ^ 

we  see  that  vf'  ,/v*  T. 

i+r  1 ' 
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The  argument  now  proceeds  as  in  corollary  3. 


(iii)  The  concavity  argument  in  (i)  shows  that  L(g  (t)) 

a-t  -1  t ac!  , , ai-  ^ 

> a(e  -1)  /q  e L(s)ds  > -(e  -1)  thus  L(t)  >“(e 

Since  ^ an.'i  both  are  increasing  g^^(t)  > t+a  thus 

ag^^(t)  ^ ^ . at+1  , ^-1,  . . II 

-(e  -1)  t(a)  > -(e  -1;  tCa).  || 


-1 

-1)  t(a). 


>-l  -1/2 


Example.  f(x)  = (r(l/2))  x"  ^ e“  , x > 0;  this  is  the  r(l/2,l)  = X^/2 
distribution  which  is  DFR  ([5]^P*  378).  The  moment  bounds  given  in 
theorem  3 and  corollary  4 for  IMEIL  distributions  apply,  as  well  as  the 
improved  moment  bounds  for  DFR  distributions  given  in  corollary  4. 
Using  our  recursive  formula  (4)  we  compute  <^  = l/2,  *^1  = ^2  ~ 

= 15/64,  c^^  = 21/32,  c^  = 315/128,  Cg  = 1485/128.  Next 
^-1  = ^0  = ''l  = ^2  = ^5/8,  vj  = 14/5,  = 15/4,  v^  = 33/7. 


Denoting  the  lower  bound  given  in  theorem  4 by  B(t)  = U(t)-  min  c.t 

0<i<6  ^ 

we  obtain:  ~ 


-1 


r 2t+i,  0 < t < 1/4 
I 2t+3/2-(8t)“^,  l/4  < t < 1 

2t+3/2-(8t^)"^,  1 < t < 15/8 


B(t)  =/ 


2t +3/2- (15/64 1^),  15/8  < t < i4/5 
2t+3/2-(2l/32t^),  14/5  < t < 15/ii 


j 2t+3/2-(3l5/l28t^),  15/4  < t < 33/7 
(^2t+3/2-(l485/l28t^),  33’l7  < t < do 


The  lower  bound  given  in  corollary  4 lidiich  we  denote  by 

B*(t)  = U(t)  - min  ctt"^,  can  be  similarly  written. 
0<i<6  ^ 
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The  table  below  gives  a few  values  of  t along  with  the  corres- 


ponding intervals  [B(t)^U(t)]  and  [B*(t),U(t)]  for  M(t). 


‘i 

(B(t),U(t)) 

(B*(t),U(t)) 

.1 1 

[1. 2,1-7] 

[1.2, 1.7] 

.5  ^ 

[2.25,2.5] 

[2.375,2.5] 

1 

[3  .375,3.5] 

[3.444,3.5] 

1.5 

[4-.  444,1+.  5] 

[4.475,4.5] 

2 

[5.471,5.5] 

[5.488,5.5] 

3 

[7.492,7.5] 

[7.497,7.5] 

[9.4976,9.5] 

[9.4990,9.5] 

5 

[11.49926,11.5] 

[11.49971,11.5] 
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